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A mixture of heavy atoms in a Mott state and light spin-1/2 fermionic atoms is studied in an 
optical lattice. Inelastic scattering processes between both atomic species excite the heavy atoms and 
renormalize the tunneling rate as well as the interaction of the light atoms. An effective Hamiltonian 
for the latter is derived that describes tunneling of single fermions, tunneling of fermionic pairs and 
an exchange of fermionic spins. Low energy states of this Hamiltonian are a Neel state for strong 
effective repulsion, dimer states for moderate interaction, and a density wave of paired fermions for 
strong effective attraction. 

PACS numbers: 03.75.Ss, 03.75.Mn, 34.50.-s, 71.10.Fd 

I. INTRODUCTION 

Recent experimental work on ultracold bosonic and fermionic atoms opened a field to study quantum 
states in optical lattices, like superfluid or Mott states [j], 0, where the interaction between the atoms 
can be controlled by a magnetic field via a Fcshbach resonance [H . Spin-dependent effects [i[ , formation 
of dimers from fermionic atoms [j| , and mixtures of two atomic species 0, 0, Q provide opportunities for 
more complex quantum states. An interesting example of the latter are dimer states. They have been 
discussed in solid-state systems @ and recently also for an ultracold Bose gas with ring exchange fiol ]. 
Frustrated spin systems were discussed for a number of atomic systems [TTj . 

The idea of this paper is to study a system in an optical lattice that undergoes a first-order phase 
transition when a model parameter is changed. If this parameter is easily accessible in the experiment 
it can be adjusted exactly at the transition point. As a consequence, the quantum system at the tran- 
sition point describes a physical situation that has a number of interesting features. First, there is the 
coexistence of the phases in form of a metastable state, leading possibly to phase separation. Second, 
new translational-invariant quantum states can appear. A system that realizes this type of physics is an 
atomic mixture which is subject to ineleastic scattering between the atomic species. It can be controlled 
by adjusting physical parameters, including the optical-lattice parameters (by choosing the frequency and 
amplitude of the Laser field) and the fermion-fermion interaction through a Feshbach resonance. This 
enables us to create ground states as well as excited states of the system, and to study their dynamics. 

In the specific case of this paper, an optical lattice is filled with heavy (bosonic or fermionic) atoms 
(HA) (e.g. 87 Rb or 40 K), one at each lattice well, which form a Mott state. The tunneling of these HA is 
neglected since the potential barriers of the optical lattice are sufficiently high. Excitations are only due 
to collisions with other atoms. For this purpose light fermionic atoms (LFA) (e.g. 6 Li) are added to the 
optical lattice. These atoms can tunnel because of their low mass and scatter by the HA. It is assumed 
that the HA experience a harmonic potential in the optical lattice, at least at low energies. Then their 
excitations are harmonic-oscillator states. During the scattering process the HA can also transfer energy 
to the LFA. Moreover, we consider spin-1/2 fermions with local (on-site) repulsion. This can be described 
by a Bose- Fermi model that is known in solid-state physics as the Holstein-Hubbard model 12]. The 
excitations of the HA are local phonons. This model is quite different from models with density-density 
interaction, where only thermal fluctuations of the HA and elastic scattering between the atomic species 
were studied [l5j . 

II. HOLSTEIN-HUBBARD MODEL 

LFA with spin a (=|, |), defined by fermionic creation and annihilation operators cj^ and c- }a in a well 
of the optical lattice with position j, tunnel between nearest-neighbor wells. The corresponding tunneling 
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Hamiltonian in tight-binding approximation reads 

H t = - J E E (CjU'^ + ^O- (!) 

<j,j'> cr=T,l 

HA, forming a Mott state, are harmonic oscillators at each well with cigcnfrcqucncy cjo and energies 
En = uiqN (N = 0, 1, ...). It is assumed that a HA in one well is excited independently of the HA in 
the other wells. This gives the Hamiltonian oJo^2jbjbj, where bj and bj are creation and annihilation 
operators for the oscillator excitations at each well, respectively. They can be considered as local phonons. 
LFA are scattered by HA due to an exchange of excitations of the HA (i.e. by creating or absorbing an 
excitation), defined by the Hamiltonian 



Hr = 



[wo&j&j + g(b] + &j)(nj T + nj|) + Un i]nil . (2) 



Here a local repulsive interaction between the LFA with strength U > has been included. The phonons 
couple to the fermionic densities and with strength g. Thus the atomic mixture is defined by the 
Hamiltonian H = H t + Hi, also known as the Holstein-Hubbard model [121 ]. 

For a given ensemble of LFA, represented by integer numbers rij a = 0, 1, the Hamiltonian Hi can be 
diagonalized with product states Y[. l-^Yj! °j) (f° r Nj = 0> 1> 2, Cj = 0, j, |, []). \Nj] 0) is an eigenstate 

of the phonon- number operator bjbj with eigenvalue JVj. A lattice product of these states is an eigenstate 
of the Hamiltonian Hi. Other eigenstates of Hi with fermions can be constructed from 

\Ny,<7)=clexpl-^(bl-bj)}\Ny,0) (<r=U) (3) 

and from a state with two fermions per well 

l^it) =cJ T cJ i e X p[-2-^(6j-6 j )]|JV j ;0). (4) 

The eigenvalue of Ylj \Nj, aj) with respect to Hi is 



E [<*> N i 



g 2 

"(«jT + n ii) 2 + Un il n il , (5) 



w 
j 

where the integers rij a — 0, 1 count the number of fermions at the well j with spin a. This expression 
represents an effective attractive fermion-fermion interaction of strength g 2 /wq which competes with 
the repulsive fermion-fermion interaction of strength U . Thus the system is given by the transformed 
Hamiltonian 

H 7 = TH t T^ + coo E b l b i + ^ E "JT^ji - — E("JT + ^l) ( 6 ) 
— — u^o — 

j j j 

with the effective coupling 7 = U — 2g 2 /loq and the unitary operator 

T = n(E E e-^ b J t -^^|7V j ;a)(7V j ; ( T|] , (7) 
j \iVj>0(7=o J T ) WT / 

where ao = 0, a-\ = ai = 1, and = 2. The states cr) contain phonons even for ATj = due to the 
exponential factors in Eqs. Q and J4|. Therefore, Nj > will be called phononic excitation, in contrast 
to = 0, which is the phononic ground state. The ground states of Hj have Nj — and are Y[ |0; It) f° r 
U < 3g 2 /uo and n(°i|0; i) + a 2 10; f)) for J7 > 3g 2 /uJo with 1 0,1 1 2 + | a 2 1 2 = 1- The corresponding energies 
are Eq = U — 4g 2 /luq and Eq = —g 2 /ujo, respectively. For strong repulsive interaction there is only one 
fermion in each well, whereas a strong attractive interaction leads to a formation of a local fermion pair 
(i.e. a bipolaron (H)])- Both interactions are balanced for U c = 3g 2 /u)q, leading to a degeneration of the 
two states. The ground-state energy Eq as a function of U has a cusp at U = U c . 

It is crucial that unpolarized spin 1/2 fermions are considered. In the case of spin-polarized fermions 
(e.g. in a magnetic trap), the only local interaction is due to Pauli's principle. This excludes multiple 
occupation of a well by LFA. Consequently, the coupling to the HA results only in a (/-dependent chemical 
potential // = —g 2 /uJo because of n? = rij in Eq. (J5]). This favors the complete filling of the optical lattice 
with one LFA per well. 
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A. Half-filled States 



Now it is assumed that the atomic system is prepared such that the number of LFA is equal to the 
number of lattice wells (i.e. 2M fermions for 2M lattice wells). This is the case of half filling and it is of 
particular interest because it allows highly degenerate quantum states. The energy of THfT^ for a half- 
filled system is calculated as a state that can have 2j singly occupied wells with energy E\ — — 2(g 2 /uj )j 
and M — j doubly occupied wells with energy E 2 = (U — Ag 2 /ui )(M — j). The total energy then is 

E a = eM - jj with e = U- 4g 2 /^o • (8) 

Here 7, the effective fermion-fermion interaction of Eq. ([6]), is an important parameter that controls two 
different phases: For 7 > the ground state has j — M (i.e. 2M singly occupied wells), whereas for 
7 < the ground state has j = (i.e. M doubly occupied and M empty wells). The ground-state energy 
is — 2(<7 2 /wo)M for 7 > and eM for 7 < 0, respectively. At 7 = the ground-state energy E does 
not depend on the number of singly occupied wells, i.e. it is highly degenerate. This point describes a 
first-order phase transition in terms of THjT^ , changing from a singly occupied lattice for 7 > to a 
mixture of doubly occupied and empty wells for 7 < 0. Moreover, both states are highly degenerate even 
for 7 7^ because the local Hamiltonian Hi does not determine how the spins, the empty, and the doubly 
occupied wells are distributed in the lattice. At U = U c a pair of singly occupied wells can be replaced 
by a pair of a doubly occupied and an empty well without costing energy. These degeneracies are lifted 
by the tunneling term H t if the more realistic Hamiltonian H — H t + Hi is considered. It is known that 
this leads to a staggered spin state (Neel state) if U U c and to a state with alternating empty and 
doubly occupied wells if U <C U c [12J. The latter is a density wave of fermion pairs. This is plausible if a 
nearest-neighbor pair of wells is considered: Two tunneling processes are possible in the case of a pair of 
singly occupied wells with opposite spin orientation but due to Pauli's principle no tunneling if the spin 
orientation is equal. Two tunneling processes are also possible in the case of a nearest-neighbor pair of 
an empty and a doubly occupied well. On the other hand, tunneling is blocked by Pauli's principle for a 
nearest-neighbor pair of doubly occupied wells. Since tunneling lowers the energy, the Neel state is close 
to the ground state for U U c and the density wave is the ground state for U <C U c . In the remainder 
of this paper the situation of U ~ U c shall be studied. 



B. Partition Function and Resolvent 



Many physical quantities of an atomic system at temperature T = 1/(3 can be calculated from the 
partition function 

Z = Tre-? H = I Tr [{z - Hy 1 ] e "^— (9) 

which is expressed here by the resolvent G{z) = (z — H)^ 1 . T is a closed contour in the complex 
plane that encloses all eigenvalues of H and Tr is the trace over all quantum states of the system. At 
low temperatures this expression is dominated by the lowest eigenvalues of H, and one can use the 
approximation Z m Zq with 



Z = ^Tr[P (z-H)- 1 P ]e-^^, (10) 

where the trace is projected by the projector Pq to quantum states that are energetically close to the 
ground state of the Hamiltonian H. In the specific case of H = H t + Hi the restriction to states without 
phononic excitations (i.e. iVj = 0) 

p ° = n(Ei°; CT i>M) 

is an obvious choice for sufficiently large ujq (i.e. sufficiently strong amplitude of the optical lattice). 
However, virtual phononic excitations are possible, excited by the Hamiltonian of er@ H in the partition 
function. 
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As an example, we consider the case without tunneling (J = 0). The projection Pq is onto half-filled 
states. The poles of the projected resolvent then read Zq = eM for a density- wave state and z = (e— j)M 
for a single-fermion state. Depending on the parameter 7, the partition function Zq is dominated by the 
lowest pole: 

-^log(Zo) = 6-76(7) , 

where 6(7) is the Heaviside function. At the transition point 7 = this expression is continuous in 7 
but its first derivative jumps with respect to 7. According to the Ehrcnfcst classification scheme of phase 
transitions, this behavior characterizes a first-order phase transition. 



III. RECURSIVE PROJECTION APPROACH 



In order to evaluate the Po-projected resolvent Gq(z) — Pq(z — H)~ 1 Pq a recursive projection approach 
can be applied [13j| . The main idea of this method is to control the mapping of states |^') = P| V I'), 
provided by the Hamiltonian H. For an initial state |\Po) a sequence of linear independent states \^k) 
(k = 1, 2, ...) is created by a successive application of H . This idea has been used in the Lanczos method 
and in the doorstep approach fl4| . The sequence of states is related to a sequence of projectors P&, acting 
on the Hilbert space of all quantum states of the system. For an initial projector Pq and its complement 
Pi = 1 — Pq this sequence is created by the Hamiltonian H . For instance, projector P2 is defined by 
applying H to Pq and projecting with Pi: Pi is the projector of the space that is created by the operator 
P\HPq, This can be repeated by using P 3 HP 2 to define P4, where P 3 = Pi — P 2 etc. It is important 
that the projectors with even index Py project onto separate subspaces. An advantage of the method is 
that for a proper choice of Pq the spectrum of the projected Hamiltonians is shifted to higher energies at 
each step. In the specific case of this paper the shift will be provided by phononic excitations. With the 
projections a sequence of projected resolvents is defined as 

Go = Pq{z-H)- 1 Pq, 



G 2 j+ 2 — P2j+2 (P2j + l(z — H)P2j+l) Pij+2 

for j = 0, 1... which are connected by the identity (more details in Ref. [l3l ]) 



G 2j = (z- P 2j HP 2j - P 2 ,PG 2 , +2 PP 2j )- i . (12) 

The iteration terminates at the j th level when P 2 j projects onto a subspace whose basis consists only 
of eigenstates of H. Such a subspace we call eigenspace. Moreover, when P 2 j projects on a subspace 
that is close to an eigenspace (i.e. a basis of this subspace has very large overlap with eigenstates), the 
truncation of the iteration is a good approximation. In terms of the Hamiltonian this means that HP 2 j 
and P 2 jP are small. This fact allows us to control the quality of the approximation. 

Here only the lowest poles are of interest because they dominate the partition function in Eq. (jXDJ) . 
Depending on the Hamiltonian H and the choice of the initial projection Pq, the lowest pole Zq of the 
projected resolvents G 2 j (z) can be increasing with j. If the values of z are restricted to those where Go (z) 
has poles but G 2 (z) has not, the sequence of projected resolvents can be truncated by approximating 
G 2 (z). The recursive projection method, given by Eq. (|12[) . is now employed for the Hamiltonian P 7 
of Eq. ©• This Hamiltonian describes tunneling of LEA in the optical lattice and interaction with HA, 
and direct interaction between the LFA. 

The initial projector Pq is given in Eq. (jlip . This implies 



P = 1 - Po 



{JVj>0} 



n 



0"; 



where restricts the summation to iVj > 1. The p)-projected Hamiltonian experiences a renormal- 
ization of the tunneling rate: 



PqH 1 Pq — Pq 



j 



5 



= (e - 7 )Af + e- 92/uJ <'P H t P + 7 ^P n jT nj i Po . (13) 

.i 

An important consequence of the coupling between HA and LFA is the renormalization of the tunneling 
rate as 

J -> T = e- g2 ' u «J , (14) 
which is the well-known polaron effect [l2j • 

A. Effective Hamiltonian 

The recursive projection will be applied to a restricted Hilbert space, where the recursion terminates 
at level j = 2. This should be sufficient for small tunneling (i.e. for J <C ojq)- Then the resolvent Go can 
be written as Gq — (z ~ H ) _1 with the effective Hamiltonian 

H = P HP + P HG 2 HP , (15) 

and G2 = [z — H%) with 

H 2 = P 2 H I P 2 + P 2 H t P 2 + P 2 HG 4 HP 2 . (16) 

The first term is 

P 2 P/P> = (e - 7)Af + u p A b i P2 + T E P 2 n jT"jl p 2 , (17) 

j j 

whereas the third term vanishes because PP2 produces only states outside of the restricted Hilbert space. 
Moreover, it is assumed that the second term vanishes because H t P 2 produces only states outside of the 
P2 -projected Hilbert space. This allows us to evaluate G2 and subsequently the effective Hamiltonian 
Hq in Eq. (fT5|) . The first part of Hq was given in Eq. [)13p. and the second part PqHG 2 HPo must be 
calculated now. The latter separates into two parts, namely a term H s that exchanges single fermions 
with opposite spin and a term H p that exchanges a fermionic pair with an empty well: 

^= t2 E E Po^Cy aG 2 cl, a ,C ia ,P EE £ ftxjj, (18) 

<jj'>T : <T' = T,i <jj'> 

describes the exchange of spins: ft-ijj' gives Q,T) < — ► (Til) f° r c ^ cr' and is diagonal for cr = cr'. 
Moreover, 

^ = t2 E E p o C ] CTC j' CT G 2C ] CT , Cj ^Po = X! ft 2jj' ( 19 ) 

<jj'> ct,<j'=T,1 <j,j'> 
describes the tunneling of two fermions as the exchange of a doubly occupied well with an empty well 

(<UT)<— » (IT,o). 

B. Eigenstates for two wells 

Now the effective Hamiltonian Ho for two wells is studied, where 

Ho = e - 7 + fttjj/ + hijj/ + /i 2 jj' + 7 (^jT«ji + "j'T n J'l) ( 20 ) 
with the single fermion tunneling term 

fc*JJ' = ~ r E ( C jV c J' CT + c \<y c i°, 



6 



Eigenstates of H on the Po-projected Hilbert space can be constructed from dimer states: 

ki'> = (IU>-|U»/V2, |d jJ ,) = (|4t,o) + |o,4t))/V2- (21) 

These are eigenstates of h\\ t y with eigenvalues Ki and 

ftxjJ'l*JJ'> = #i|*jJ'> , h 1<jt y\d^)=0 , (22) 
and also eigenstates of hi jj' with eigenvalues K% and 

^JJ'IdjJ') = foldjJ') , hajj^jj') = . (23) 
The eigenvalues, which have been plotted in Fig. 1, are 

K 1 = 2T2 yi- ( 2 *>° 2 > m , ^^vy- 2 ^" . (24) 

z — ' to! 2 — 2e + 27 — ajo?7i to! z — 2e — worn 

m>l m>l 

hkjj 1 (k = 1,2) do not connect different dimer components: 

ftfcJJ'UT,*) = VjJ'kat) =^.j,j'|0,a) = /*jj/|M) =0 (<r =i, T; fc = l,2). (25) 

Moreover, ht,yy switches between dimers with rate 2t: 

h t,i,y\s},y) = 2r|djj/) , Vjj'Mjj') = 2r l s j,j'> ■ 

Since always i^i < (cf. Fig. 1), the formation of dimers |sy) is favored. However, when 7 < the 
disadvantage of the dimer \dyy) can be compensated by the local interaction term in H$ which is diagonal 
and counts the number of doubly occupied wells. For the state a\syy) + b\dyy) we get the relation 

H (a\s 3S ) + &|d,j,» = (e - 7 + K x + 2rb/a)a\syy) + (e + K 2 + 2ra/b)b\dyy) . (26) 

This state is eigenstate of hyy for 7 = K 1 — K2 + 2r(b 2 — a 2 )/ab with eigenvalues e — 7 + K 1 + 2rb/a. 
For a given 7 we can evaluate o and b and then determine the related eigenvalues. In particular, for 
7 = K\ — K2 there are solutions a = ±b with corresponding eigenvalues e + K% ± 2r. 

Ki, K2 still depend on the parameter z which has to be determined as a pole of the resolvent Gq(z): 

z = 6 + K 2 (z)±2t , (27) 

where the lowest pole is the most relevant contribution. Other solutions of this equation refer to excited 
states that are not relevant for the discussion here. 

C. Eigenstates for the lattice 

On the lattice with 2M wells the effective Hamiltonian Ho can be constructed from the effective 
Hamiltonian in Eq. (|2H)) by a summation over nearest neighbors. In order to study simple eigenstates of 
the lattice system first, the tunneling of single fermions is excluded from the discussion: 

H =H - £ h t ,yy . (28) 
<jj'> 

Eigenstates of H' are product dimer states 



with eigenvalues 



i*i> = n m . i*2> = n (29) 
<jj'> <jj'> 



Ar = (e- 7 )M + 2Mr 2 £ -1 ^'ff , (30) 

' to! z — e — 7)iW — 
m>l \ u u 
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a 



FIG. 1: Behavior of K\ (full curve), K2 (dashed curve) and —r (dotted curve) vs. the strength of the fermion- 
phonon coupling a — g/uio. The other parameters in K\ and K2 are wo = J= l,7 = 0, and 2 is given by the 
lowest pole. All energies are measured in units of loq. 



and 

„2 /,.,2\m 



A 2 = eM + 2Mr 2 V -1 ±^IMT_ (31) 

m>l 

The alternating dimer state 

\*ad) = J] ( a JJ'|ajJ'> + fcjJ'MjJ'» 
<j,j' >ec 

is also eigenstate of H' with eigenvalue 

\ AD = (e-j/2)M + Mr 2 V -1 ^fe^ + Mr 2 V -1 ■ (32) 

v " ; ^ m! z- e -7/2 )M- w m ^ m! z — (e — j/2)M — uj a m v ; 

m>l m>l 

The coefficients aj.j' and fry are either or 1, alternating on the lattice. C covers the lattice with dimers, 
separated by nearest-neigbor bonds (cf. Fig. 2). 

Next, the energy z must be determined by the lowest poles of the resolvent Go- This gives, up to small 
corrections, for Ai the value z = (e — j)M such that 



r 2 



Ai/M = e - 7 - 2— F(. 9 2 /^o) with F ( x ) = ' (33) 

LOq m >l m ' m 

For A2 the value z — eM gives 

A 2 /M = e - 2—Fi-g 2 /^) , (34) 
and for Xad the value z = (e — j/2)M leads to 

2 2 

X AD /M = e - 7/2 - -F( 3 2 /^ 2 ) - ^F(- 3 2 /^ 2 ) . (35) 
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4^4 1=t 

FIG. 2: The alternating dimer state of Eq. (|21|l on a square lattice, where the elements represent | J,, f) — | f, J.) 
and | |t, 0) + |0, ||), respectively. 

Dimer states in a bipartite lattice are not expected to be ground states. For instance, the Neel state gives 
for z — (e — j)M 

2 

(* N \H \* N )/M=(e- 1 ) + {H> N \ V h 1 ^ N )/M=(e-j)-^F(g 2 /^) , 

<jj'> 

where Z£ is the number of nearest-neighbor wells. The corresponding matrix element for a density-wave 
state of fermion pairs is e, which is very close to the eigenvalue A2/M of the paired fermion dimers {^2)- 
The single-fermion tunneling has not been included, since the renormalized tunneling rate r is very 
small for sufficiently large coupling g/ujo- Therefore, its effect on the eigenvalues will be neglegible, 
provided the system is not at a degenerate point. Although the dimer states are not ground states, the 
system can stay for very long times in them, in terms of a typical experiment. 




IV. DISCUSSION AND CONCLUSION 



An effective Hamiltonian Hq has been derived within a truncated recursive projection approach for a 
half-filled optical lattice. The truncation should be valid under the assumption that the tunneling rate 
J is small in comparison with the local potential of the HA (i.e. for J <C ujo). The properties of the 
effective Hamiltonian in Eq. (p~5|) can be summarized as follows: (I) Ho acts in a four-dimensional space 
at each well of the optical lattice, spanned by {|0), | j), | |), | It)}- (H) Ho describes one and two fermion 
processes between nearest-neighbor wells. Phonons (i.e. excitations of the HA) are created by the unitary 
operator T of Eq. ([7]) and appear in the eigenvalues of the effective Hamiltonian Ho as excitation of m 
phonons. In the case of a system with two wells the eigenstates of H are superpositions of dimer states, 
where either the fermions sit in separate wells or as a pair in only one well. The contribution of each 
of the two dimers to the eigenstate can be controlled by the parameter 7 = U — 2g 2 /ujo- This result 
was extended to an optical lattice with 2M wells. If single-fermion tunneling is neglected in Ho (this is 
a good approximation for sufficiently strong fermion-phonon coupling g/cu), dimer states are eigenstates 
of the effective Hamiltonian. Although the dimer states are not ground states of the system, they can 
be prepared by using a modulated optical lattice, where the tunneling barrier is higher between the 
dimers. This is possible because we have control of the optical lattice (Laser frequency and amplitude) 
and of the fermion-fermion interaction U through the Feshbach resonance. After the preparation of the 
dimer state, the modulation can be turned off without destroying the dimers, since it is (almost) an 
eigenstate. The main effect of the (weak) single-fermion tunneling term is to flip singlett dimers into 
dimers of paired fermions (and vice versa) with a rate 2r. Therefore, single-fermion tunneling favors the 
formation of neighboring pairs of different dimers. Three different dimer states were studied: states with 
singlett dimers, states with dimers of paired fermions and states with both dimers, placed in alternating 
order in the lattice. Their corresponding energies (given by their eigenvalues with respect to the effective 
Hamiltonian) are different. However, there is a point of degeneracy for 

7c = 2-[F(-. 9 2 K 2 )"^(5 2 K 2 )]] . 



as plotted in Fig. 3. 
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FIG. 3: Eigenvalues of H' of Eq. (|28[) for dimer states: Xi/M (full curve) for singletts, X2/M (dashed curve) 
for paired fermion dimers, Xad /M (dotted curve) for an alternating dimer state vs. the strength of the fermion- 
phonon coupling a = g/u>o and fixed J = 1, U = 2. For comparison, the diagonal matrix element with respect to 
the Neel state is shown for zl — 4 by the thin curve. All energies are measured in units of ljq. 



The interaction between LFA and HA in an optical lattice can be understood in a simple physical 
picture. Already in the absence of tunneling, the local repulsive interaction U of the LFA is affected by 
the coupling to the HA. It results in a reduction of U to 7 = U — 2g 2 /uq. The reason is that the HA 
absorb interaction energy from the LFA by going into an excited state (cf. Eq. lO). I n the presence 
of a tunneling rate J of the LFA also kinetic energy is taken away by the HA, leading to a reduction of 
the tunneling rate to r = Jexp(— g /oJq), The renormalization of the tunneling rate and the interaction 
strength depends on the parameter g/ujQ. This means that the effect of the HA on the LFA is the 
stronger the easier it is to excite the HA with energy u>q. The latter can be controlled in an optical 
lattice by tuning the shape of the potential at the bottom of the wells. Since excitations of energy loq are 
essential, a deviation of the potential from the assumed harmonic form at higher energies does not affect 
the qualitative picture. The fact that 7 can be negative (attractive interaction) implies a transition from 
a repulsive to an attractive gas of LFA. At half filling the ground states of the LFA in an unfrustrated 
lattice is a Neel state in the case of strong repulsion and a density wave of paired fermions in the case 
of strong attraction. For 7 = there is no effective interaction between the LFA. However, this does 
not mean that there is a free Fermi gas because the scattering with the HA still affects the tunneling 
processes. Our calculation, based on a small tunneling rate, shows that the LFA form dimer states for 
7«0, where either two fermions on nearest neighbor wells with opposite spin orientation or a locally 
bound pair of two fermions and an empty well are formed. 

The dimer states, as well as the Neel and the density- wave state, should be observable in a Stern- 
Gerlach experiment Q. For this purpose, a mixture of two atomic species, constituted by the same 
number of LFA (e.g. 6 Li) and HA (e.g. 87 Rb or 40 K) is brought into an optical lattice, where half filling 
is maintained at low temperature. Then the interaction parameter 7 = U — 2g 2 /ujq is adjusted, either by 
tuning U through a magnetic field near a Feshbach resonance, or by tuning the excitation energy luq of 
the HA (i.e. the tightness of the wells in the optical lattice). Moreover, the tunneling rate J is adjusted 
by the amplitude of the optical lattice. After the creation of the quantum state, the optical lattice is 
turned off to allow the gas to expand and to measure the spin projections in a Stern-Gerlach experiment. 
The Neel state, consisting of individual fermions with alternating spin orientation, would separate in the 
inhomogeneous magnetic field into two parts. In contrast, the density wave of paired fermions and the 
dimer states would not separate so easily. 

In conclusion, a mixture of heavy atoms and light fermionic atoms has been investigated. Inelastic 
scattering processes generate an attraction between the light atoms and a strong reduction of the tunneling 
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rate. Depending on the total interaction, there is the possibility of a Neel state, a density wave of paired 
fermions, and different dimer states in an optical lattice. 
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